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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2939 

OPTIMUM  CONTROLLERS  FOR  LINEAR  CLOSED-LOOP  SYSTEMS 

By  Aaron  S.  BoksenLom,  David  Novik 
and  Herbert  Heppler 


SUMMARY 

An  analysis  for  optimum  controllers  of  general,  linear,  time- 
invariant  multiloop  systems  is  presented.  Optimization  is  based  on 
minimizing  mean-square  or  integral- square  errors  for  either  stationary 
statistical  or  transient  inputs,  with  limits  and  constraints  of  meain- 
or  integral- square  form.  General  representations  of  controller  char- 
acteristics, stable  process  characteristics,  and  error  relations  are 
used.  A method  is  shown  of  assuring  stability  of  the  multiloop  system 
during  the  optimization  process  and  casting  the  multiloop  controlled 
system  into  an  equivalent  open  loop  so  that  the  methods  of  optimum  fil- 
ter theory  can  be  used.  General  solutions  are  obtained  for  four  special 
cases.  Two  examples  for  speed  control  of  a turbojet  engine  illustrate 
the  methods  developed  in  this  report.  The  data  for  a controlled  turbo- 
jet engine  were  in  substantial  agreement  with  the  theoretical  results. 


INTRODUCTION 

One  of  the  fundamental  functions  of  controllers  is  the  reduction 
of  certain  variables,  called  errors,  to  small  values.  This  must  be 
accomplished  with  the  dynamical  nature  of  the  system,  the  difficulty 
of  accurate  measurement,  the  power  requirements  to  manipulate  valves, 
and  so  forth,  the  random  effects  at  the  inputs,  and  the  general 
effects  of  noise  taken  into  consideration.  At  various  stages  in  the 
design  of  an  over-all  controlled  system,  the  question  arises  as  to 
what  computations  the  controller  should  make  to  minimize  the  errors. 

In  this  report,  the  problem  of  optimizing  this  computational  aspect 
of  the  controller  is  analyzed.  Linear  systems  are  assiomed  throughout. 
Optimization  is  considered  under  either  transient  or  statistical  inputs 
and  disturbances  and  is  based  on  the  minimization  of  mean-square  errors 
or  integral- square  errors.  Constraints  and  limits  are  included  as  mean- 
square  or  integral-square  values. 
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Most  investigations  in  this  field  optimize  the  controlled  system 
by  adjusting  certain  parameters  of  an  otherwise  fixed  system.  In  this 
report  the  more  general  approach  of  optimizing  the  system  by  adjusting 
the  entire  frequency-response  characteristics  of  the  controllers  is 
used.  This  general  approach  was  first  developed  in  reference  1 for 
the  design  of  filters,  and  a corresponding  development  is  presented 
herein  for  a closed-multiloop  controlled  system.  New  problems  of 
structural  stability  and  physical  realizability  arise  in  the  closed- 
loop  cases.  In  addition,  necessary  constraints  and  limits  are  included 
in  this  analysis. 

The  problem  of  optimizing  the  computational  aspect  of  the  control- 
ler is  first  reduced  to  a standard  form  by  generalizing  the  frequency- 
response  characteristics  of  both  the  process  to  be  controlled  and  the 
controller,  the  conditions  of  stability,  and  the  specifications  on 
errors,  constraints,  and  limits.  The  general  solution  for  the  optimum 
controllers  giving  an  absolute  minimization  of  errors  is  derived.  In 
addition,  the  ejqjressions  for  the  minimum  errors,  and  the  additional 
errors  suffered  when  nonoptimum  pontrollers  are  used,  are  shown. 

It  is  expected  that  the  resixlts  of  this  analysis,  conducted  at 
the  NACA  Lewis  laboratory,  can  be  used  as  a basis  for  controller  design, 
as  a standard  under  which  controllers  can  be  evaluated,  or  as  a possible 
basis  of  specifications  on  controller  dynamics.  Applications  of  the 
results  of  this  analysis  to  several  examples  of  controlling  gas  turbine 
engines  are  shown. 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 
a engine  or  fuel  servo  dead  time 

b,c  constants  for  describing  spectral  densities  of  inputs 

C general  linear  operators  representing  controller 

E stable  linear  operators  representing  engine  or  process 

e errors  to  be  minimized  by  control  and  constrained  variables 

e*  errors  resulting  from  nonoptimum  controller 

F,G  stable  linear  operators  representing  parts  of  controller 

H general  linear  operators  giving  response  of  errors 
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K loop  gain 

R arbitrary  stable  linear  operators 

t time 

w inputs  and  disturbances  affecting  errors 

X manipulated  variables 

y Inputs  to  controller 

z inputs  and  disturbances  affecting  y's 

6F  arbitrary  stable  linear  operators 

X Lagrangian  multiplier 

engine  time  constants 
integral  time  constant 
(D  frequency 

Subscripts: 

J^n^r^s^v  summation  indices 
opt  optimum 

Superscripts: 

complex  conjugate 

expected  value  or  integral 

Correlation  function  notation: 

For  statistical  inputs, 

1 -icDu 


(WjZv)  = (zv^i)  = 


v"j‘ 


2jr 


Zy(t)w. (t  + u) 


du 


V,yU  = -ao 


For  transient  inputs, 
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n» 


w 


i(t)  e 


-ioit 


dt 


K„JLx3 


Zv(^) 


+io)u  du 


Matrix  notation: 


E., 

Jk 


indicates  element  of  row,  column  of  matrix  E 


ANALYSIS 

Scope  of  Analysis 

The  scope  of  the  analysis  is  summarized  as  follows: 

1.  Time -invariant  linear  systems  are  assumed  throughout. 

2.  The  process  to  he  controlled  is  inherently  stable,  with  one 
variable  to  be  manipulated  by  the  controller. 

3.  Cptimization  is  to  be  realized  for  either  transient  or  station- 
ary random  inputs. 

4.  Optimization  is  based  on  minimizing  mean- square  or  integral- 
square  errors.  The  constraints  and  limits  of  the  system  are  included 
by  constraining  the  mean-square  or  integral- square  values  of  the  vari- 
ables involved. 

5.  Complete  freedom  in  computation  is  assumed  for  the  controller. 

The  same  formal  equations  in  the  time  domain  are  obtained  for  the 
nonstationary,  time-varying  cases  as  for  the  stationary,  time -invariant 
cases  (ref.  2);  but,  even  for  the  open- loop  filter  problem,  the  equa- 
tions for  the  time-varying  cases  can  be  solved  only  by  numerical  means. 
Extension  of  the  analysis  of  this  report  to  the  time-varying  cases  would 
require  such  techniques  as  those  in  reference  2.  The  linear  assumptions 
do  not  completely  preclude  application  of  the  methods  of  this  report  to 
nonlinear  systems.  Many  nonlinearities  can  be  handled  so  as  to  allow 
linear  techniques  without  omitting  the  essential  nature  of  the  nonline- 
arity (refs.  3 and  4). 

When  the  process  to  be  controlled  is  inherently  unstable,  the 
techniques  of  preserving  stability  of  the  closed-loop  system  are  more 
involved  and  these  rather  uncommon  cases  are  not  treated  herein.  The 
extension  of  this  analysis  to  cases  having  several  manipulated  variables 
follows  in  a straightforward  manner,  but  the  final  equations  to  be 
solved  are  more  complex.  The  fundamental  development  of  these  cases 
is  shown  in  the  appendix. 
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The  optimization  problem  for  transient  inputs  is  the  same  as  that 
for  statistical  inputs  (ref.  5)  if  mean  squares  are  replaced  by  integral 
squares  and  if  spectral  densities  are  replaced  by  products  of  Fourier 
transforms,  as  shown  in  the  list  of  symbols.  The  analysis  which  follows 
was  guided  by  the  desire  to  keep  the  systems  studied  and  the  problems 
resolved  as  general  as  possible  so  that  wide  applications  may  be  made 
of  the  results.  In  order  to  utilize  available  open-loop  theory,  it  was 
necessary  to  determine  an  open  loop  equivalent  to  the  actual  closed 
loops  which  would  still  assure  structural  stability  of  the  actual  sys- 
tem during  the  optimization  procedure.  Techniques  of  accomplishing 
this  end  were  found  for  those  cases  of  stable  process  characteristics. 


Characterization  of  Process 


The  general  form  of  the  systems  considered  is  shown  in  figure  1. 

The  time -invariant  linear  process  to  be  controlled  has  certain  inputs 
X (in  this  case  only  one,  such  as  a fuel-valve  position  signal)  which 
can  be  and  are  chosen  to  be  manipulated  by  the  controller.  There  are 
a certain  number  of  outputs  y (such  as  measured  engine  speed,  meas- 
ured engine  speed  error,  throttle  position,  altitude  pressure,  etc.) 
which  can  be  and  are  chosen  to  be  used  by  the  controller.  Each  output 
is  affected  by  transient  and  statistical  disturbances  and  inputs  z of 
altitude,  air  speed,  throttle  position,  and  so  forth,  and  by  the  manipu- 
lated variable  tlirou^  a linear  operator  E.  In  general,  E would  be  a 
rectangular  array  of  operators,  but  for  the  case  of  one  manipulated 
variable  E is  a column  of  operators.  The  general  form  of  the  process 
to  be  controlled  can  be  written  operationally  as 


This  form  follows  directly  from  one  method  of  obtaining  these  char- 
acteristics. If  X is  held  constant,  measured  y gives  the  nature  of 
z.  If  X is  varied  sinusoidally,  the  harmonic  analysis  over  a large 
number  of  cycles  of  measured  y filters  out  any  effects  of  z because 
the  harmonic  content  at  any  one  frequency  of  either  a transient  or  sta- 
tistical z is  relatively  small.  This  can  be  seen  from  the  following 
analysis: 


If 


then 


x(t)  = Re  (e^^)  t >0 


- 0 


t <0 


y(t) 


z(t)  + Re 


as  t-*’  * 
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internal  disturbances  of  such  structvires  can  be  considered  as  applied 
to  the  input  or  output  as  small  disturbances,  the  proper  discussion  of 
stability  should  deal  only  with  stable  structures. 

The  controller  operational  equation  should  therefore  be  written 


X 


yj  + 0 . X 


(3) 


where  Fj  and  G represent  stable  structures.  The  only  new  modes  of 

oscillation  in  the  system  in  addition  to  the  stable  modes  generated 
by  E,  F,  and  G are  those  represented  by  the  response 


1 - G - 


This  criterion  of  stability  can  be  seen  by  deriving 


the  responses  of  all  variables  in  the  closed-loop  system  to  all  possible 
disturbances.  Each  response  will  have  the  preceding  factor  multiplied 
by  E,  F,  or  G,  or  by  the  sum  of  the  products  of  these  operators.  It 
can  also  be  shown  that  no  loss  of  generality  is  caused  by  allowing 


G = 


(4) 


in  which  case  stability  of  the  entire  closed-loop  system  is  assured. 
The  number  of  functions  F still  to  be  determined  is  the  same  as  the 
original  niunber  of  unknown  functions  C of  equation  (2). 


There  are  alternative  ways  of  obtaining  the  same  over-all  opera- 
tion ejq)ressed  by  equations  (2)  or  (3),  some  of  which  may  be  unstable. 
One  general  method,  always  stable,  is  shown  by  figure  2,  in  which  the 

F.’s  are  arbitrary  stable  structures. 

J 

Combining  equations  (l),  (3),  and  (4)  gives 


X 


(5) 


The  original  closed-loop  system  of  figure  1 is  now  equivalent  to  the 
open-loop  system  of  figure  3,  and  this  open-loop  system  is  now  compar- 
able with  the  systems  of  reference  1. 
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An  expression  for  Cj  of  equation  (2)  is 

F. 


z 


1 ^ N EjFj 


(5a) 


Specifications  on  Controller 

The  errors  to  he  minimized  hy  the  controller  are  characterized  in 
the  same  way  as  the  y’s  as  they  have  a controlled  and  an  uncontrolled 
part.  The  specifications  are  written 


^ • • 

e^"^  = a minimum 


J 


(6) 


where 


= ''j  + Hj  • X (7) 

and  H is  any  general  linear  operator. 

The  X multiplier  technique  is  used  to  allow  additional  degrees 
of  freedom  whenever  the  optimum  values  of  the  errors  are  not  independent 
or  when  constraints  of  mean-square  form  are  to  he  imposed.  The  X's 
may  then  appear  as  parameters  of  the  controller  and  are  adjusted  to  give 
a compromise  among  the  various  dependent  errors  or  to  set  the  mean- 
square  values  of  the  constrained  variables. 


Minimization  of  Errors 

From  equations  (?)  and  (s)^  the  expressions  for  the  quantities  to 
he  minimized  or  constrained  are 


e 


J 


z 


n 


and 
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and 


where 


and 


T 


2rtk 

ID 


(k  is  an  integer) 


L -*•* 


The  first  term  on  the  right  side  of  the  preceding  equation 
approaches  zero  as  k-^®  for  either  a transient  or  statistical  z, 
and 

L+T 

g,-io)t  for  L;  k-+® 


Characterization  of  Controllers 

The  controller  to  he  designed  is  to  operate  on  the  y's  to  give 
X according  to  the  operational  equation 


X 


(2) 


As  the  optimization  of  only  the  computational  aspect  of  the  controller 
will  he  considered^  the  y’s  end  the  x are  outputs  and  input,  respec- 
tively, which  isolate  the  engine  or  process  to  he  controlled  end  which 
can  he  freely  used  in  computation. 


The  C’s  of  equation  (2)  are  arbitrary  linear  operators  which 
need  only  to  keep  the  closed-loop  system  stable.  But,  in  order  to 
assure  stability,  the  responses  of  all  variables  in  the  system  to  all 
possible  disturbances  must  be  considered.  These  include  variables  and 
disturbances  internal  to  the  C and  E boxes.  Because  arbitrary 
linear  operators  can  be  built  up  from  stable  structures  and  the  small 
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^^|Hj|2f^^(z^z^)  + (wjVj) 
V n - 


(8) 


If  the  symbol  F is  used  for  the  optimum  system  resulting  in 
error  e and  the  symbol  F + 5F  is  used  for  any  other  system  resulting 
in  error  e*,  then 


+ 


H 


^n^vJ 


n 


+ 


The  first  two  terms  of  the  right-hand  side  of  the  preceding  equa- 
tion are  equal  because  the  integrands  are  merely  conjugate.  The  third 
term  is  nonnegative  as  it  is  equal  to  the  mean- square  error  resulting 
when  F is  6F  and  v ~ 0.  Thus  the  necessary  and  sufficient  condi- 
tion for  Ein  absolute  minimum  is  that 


j 


1 da)  ^ 5Fy 

Hj(z^wj)  + ^IHjI^F^(z^z^) 

- n _ 

(9) 


eind  then 
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opt 


fZ  l"j  Z Z ‘ *0 

(10) 


Because  6F^  is  any  arbitrary  but  structurally  stable  system, 
equation  (9)  reduces  to 

^ XjHj  (z^Wj)  + j ^ \j  |Ej  j 2 j ^ F^(z^z^)  = ^ 

j \ J y n 


(11) 


for  all  values  of  v for  which  there  are  a corresponding  Fy  and  Cy, 
and  where  represents  any  arbitrary  stable  structure. 


Cases  in  Which  z's  Are  Independent 

A general  solution  of  equation  (ll)  has_been  found  only  for  the 
case  in  which  the  z's  are  independent  ((z^Zy)  is  a diagonal  matrix). 

In  this  case  the  general  solution  is 


Fy  = 


Vl 

MZ^ 


(12) 


where 


= _j 

■^~v  

MZ^ 

X b Ibl  ^ 

j 

(Zyly)  = iZyl^ 
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The  operators  M and  represent  minimum  phase  structures^  and 

A^  = -^,1  + ■^;2  ■where  A_^^l  and  Ay.^2  represent  stable  structures. 

This  solution  follows  the  general  pattern  of  reference  1 in  the  fre- 
quency domain.  The  determination  of  M and  Zy  is  the  factorization 

problem  of  reference  1.  The  f\mction  A^^p  can  be  obtained  formally 

by  the  following  equation: 


From  equation  (lO),  the  difference  between  nonoptimum  and  optimum 
errors  is 


dio|M5F^(^ 


(14) 


EXAMPLES 
Example  1 

The  first  example  taken  from  a problem  in  the  control  of  a turbo- 
jet engine  is  illustrated  in  figure  4.  The  engine-speed  error  is  fed 
into  the  controller  and  a fuel-valve  position  signal  is  varied  by  the 
controller.  Linearized  engine  characteristics  are  represented  by  a 
lag  (ref.  6)^  and  the  fuel  servo  is  considered  as  having  only  dead  time 
(the  output  reproduces  the  input  a seconds  later). 

Then 

6i  = yi  = engine  speed  error 
= Zp  = negative  of  speed  setting 
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,-aitb 

^1='^1  = i +xl£0 


X = signal  to  fuel  servo 


From  equation  (l2),  since  |M|^ 


then 


Letting 


then 


|l  + TiO)  p ^ 


1/2 


M = 


1 + xiu) 


Ai  = 


(z^Zi)  = ^ 


|h  + icjo|‘ 


Zt  = 

h + 10) 


A 


■1 


= 


aio) 

e c 

b + io) 


*1,1  = s I « 

lo  ^ 


X 1/2  c 
b + iu 


u=- 


and 


rv 


-ia)t  l/2  -b(t+a) 

6 dt  Xj^  ' c e ^ ' 


Uo 

“ c C 


X,l/2  e-* 


b + io) 


-ab 


Fq_  = -e  (l  + tIo)) 
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From  equation  (5a) 


From  equation  (l3) 


- + Tio)) 

2^  _ g -atg-aio) 


-2ab^ 


For  the  controller  action  derived,  the  responses  of  speed  and  x to 
speed  setting  are 


Speed  = • (speed  setting) 


and 


X = e"^^(l  + Tun)  • (speed  setting) 

The  order  property  and  thus  the  derivative  action  of  the  preceding  equa- 
tion indicates  the  large  fuel  flow  and  temperature  variations  encoun- 
tered with  the  derived  controller  because  no  limit  or  constraint  was 
placed  on  these  variables. 


Example  2 

The  second  example  in  which  the  same  control  problem  as  in  exam- 
ple 1 is  analyzed  with  a constraint  on  the  mean-  or  integral -square 
turbine  temperature  variation  and  in  which  the  servo  lag  has  been 
omitted  is  illustrated  in  figure  5.  For  this  example, 

ej_  = y^  = engine  speed  error 
Wf  = Zi  = negative  of  speed  set 
eg  = turbine  temperature 
wg  = 0 


i i 1 + Tun 


„ 1 + aini 

^2  - 1 + TiCD 


X = signal  to  fuel  servo 
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From  equation  (l2),  since  |Mp  = 


^2  I ^ cric)  ^ 

Lt 


M = 


|l  + Tio)  1^ 
(^2^  + X2)  ^ t ^2  ^ 


1 + xiyii 


then 


Ai  = 





Letting 


|b  + id)  1 2 


then 


=F 


+ id> 


^ (h  + i(Jo)  (X^  + ^2^^^^  - X2^/^cri£nJ 


X^c 


Xj^c 


(h  + id))  (Xi  + + X2^'^^crh^ 


and 


^1  - 


X]^(l  + Tid)) 

pXi  + X2)V^  + X2^/^cjhJj|(Xi  + X2)’^'^^  + X2^/^aidJ 
From  equation  (5a) 

-Xi(l  + Tico) 

|^(Xi  + Xg)^/^  + ^^2^/^crhJQxi  + X2)^/^  + - Xi 

From  equation  (l3) 


Cl  = 
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[xgd  ^ 

h 0^b^)  + 2(X^  + X2)^/^X2^/^crb 

(X^  + X2)V2  + X2V2ab  ^ 

(15) 


The  responses  of  speed  and  x to  speed  setting  are 


Speed  = 


X = - Ft  • (speed  setting) 


(speed  setting) 


In  all  the  preceding  equations,  only  the  ratio  of  the  two  X*s  is 
effective.  This  ratio  can  now  he  set  by  evaluating  the  individual 
errors  that  make  up  equation  (l5).  Choosing  the  ratio  of  the  X's 
involves  a compromise  between  temperature  and  speed  errors,  as  these 
quantities  vary  oppositely  to  the  ratio  of  X’s. 


EXPERIMENTAL  RESULTS 

An  axial-flow  turbojet  engine  was  operated  on  a sea-level  static 
test  stand  to  obtain  some  verification  of  the  analysis.  Figure  4 of 
example  1 describes  the  system.  The  frequency  response  of  engine  speed 
to  fuel  pressure  obtained  by  harmonic  analysis  of  transient  data  is 
shown  in  figure  6.  An  approximation  to  this  data  gave  T = 1.6  seconds 
and  a = 0.163  second. 

A proportional-plus-integral  contro3J.er  of  the  form 
C = -K  ^1  + was  used  over  a range  of  K’s  (loop  gain)  from 

3.5  to  11  and  a range  of  T2^'s  from  1 to  2.5  seconds.  The  integral- 

square  percentage  error  was  obtained  for  each  controller  setting  from 
the  response  to  a step  in  speed  setting,  and  the  results  are  shown  in 
figure  7. 

The  smallest  error  (of  magnitude  0.41  sec)  was  obtained  at  K = 9.3 
and  =-1.6  seconds. 

The  results  of  example  1 with  b = 0 (step  input)  gave 


Copt  - 


+ Tunj 


1-6 


-aitjo 
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1 

With  C, 


1 + 


aW  opt 


I (i  ^ 4)(^  - ¥)■ 


The  small  deriv- 


ative action  indicated  was  not  used.  For  t = 1.6  seconds 

and  a = 0.163  second;  the  preceding  equation  indicates 


K 


opt 


1.6 

0.163 


= 9.82 


T-  , = T = 1.6  seconds 
l,opt 


which  is  in  substantial  agreement  with  the  data  of  figure  7. 


SUMMARY  OF  RESULTS 

An  analysis  was  developed  for  the  optimum  controllers  of  general 
linear  closed-multiloop  systems  under  either  stationary  statistical  or 
transient  inputs  to  minimize  mean-square  or  integral- square  errors. 
General  solutions  were  obtained  for  the  case  of  independent  inputs  and 
one  manipulated  variable.  Several  examples  of  this  case  were  shown 
for  the  speed  control  of  a turbojet  engine.  Experimental  data  from  a 
controlled  axial-flow  turbojet  engine  were  in  substantial  agreement 
with  the  theoretical  results. 

For  the  general  multiloop  systems  of  any  number  of  manipulated 
variables;  general  solutions,  which  are  shown  in  the  appendix,  were 
obtained  for  four  special  cases. 

Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland;  Ohio,  Febmary  16,  1953 
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APPENDIX  - GENERAL  MULTILOOP  SYSTEMS 

The  scope  of  this  analysis  is  the  same  as  for  the  preceding  analy- 
sis except  that  it  includes  any  number  of  manipulated  variables. 

Characterization  of  system.  - The  general  form  of  the  system  to  be 
considered  is  shown  in  figure  8.  The  general  form  of  the  process  to  be 
controlled  is  written  as  follows  in  matrix  notation: 

y = z + E • X (l6) 

The  controller  to  be  designed  is  to  operate  on  the  y*s  in  order  to 
give  the  x's  according  to  the  operational  equation 

X - C • y (17) 

or 

X = P . y + G • X (18) 

where  every  element  of  the  F and  G matrices  represents  stable 
structures.  The  new  modes  of  oscillation  generated  by  the  closed  loops 
are  represented  by  the  responses 

[l  - G - FE]'^ 

No  loss  in  generality  is  caused  by  allowing 

G = - FE  (19) 

and  stability  is  now  assured.  There  are  alternate  ways  of  building  the 
same  over-all  operation  expressed  by  equation  (l7)  or  (18),  some  of 
which  may  be  unstable.  One  definitely  stable  method  is  shown  in 
figure  9. 

Combining  equations  (l6);  (18),  and  (l9)  gives 

X - F • z (20) 

which  represents  an  equivalent  open-loop  system.  An  expression  for  the 
C of  equation  (l7)  is 

C = (l  + FE)"^F  (20a) 

Specifications  on  controller.  - The  errors  and  constraints  are 
characterized  in  the  same  way  as  y.  The  specifications  are  thus 
written 
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= a minim-um 


(21) 


vhere 


e = w + H • X (22) 

and  the  elements  of  H are  general  linear  operators. 

Minimization  of  errors.  - From  equations  (22)  and  (20),  the  expres- 
sion for  the  quantities  to  he  minimized  or  constrained  is 

e = w + HF  • z 


and 


n 


dn> 


^ (HF)jy(zyWj)  + ^ ^ 


(BF)jy(HF)jj^(z^zJ  + (wjWj) 


V n 


If  the  symbol  F is  used  for  the  optimum  system  resulting  in 
error  e and  the  symbol  F + 5F  is  used  for  any  other  system  resulting 
in  error  e*^  then 


4 


X - (X  = X X X jn("n^v) 

j J ^pt  J J.„  ^ L n J 


OO  ■ 

X do)  X (H&F)jy  (z^Wj)  + ^ 

J J-00  V L n 


(HF)jn(ZnZv) 


2 r 22  ja<  Vb> 


J V n 

^-OD 


2789 


2789 
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The  first  two  terms  of  the  right  side  of  the  preceding  equation 
are  equal  since  the  integrands  are  merely  conjugate.  The  third  term 
is  nonnegative  as  it  is  equal  to  the  mean- square  error  resulting  when 
F is  6F  and  w = 0.  Thus  the  necessary  and  sufficient  condition 
for  an  absolute  minimum  is  that 


^ Xj  do)  ^ (H6F)j^  |^(z^wj)  + ^ 


(SF)jn(2n^v) 


= 0 (24) 


and  then 


j \ J /opt  J V n 


(H5F)j^(H5F)j^/z^z^)^0 


(25) 


Because  each  element  of  the  6F  matrix  is  any  arbitrary  but  structur- 
ally stable  system,  equation  (24)  reduces  to 


= Rrv  (2^) 


for  all  values  of  r and  v running  over  the  rows  and  columns  of  F, 
and  where  every  element  of  R represents  any  arbitrary  stable  structure. 

Equation  (26)  in  matrix  symbols  becomes 

I'XHFZ  + I'XW  = R (27) 

where  H’  indicates  the  transpose  of  H and  the  matrices  Z,  W,  and 
X are 

2jv  = 

X=  diagonal  (Xj^,  Xg;  Xj,  . . .) 

The  general,  solution  of  equation  (27)  or  (26)  has  not  been  found.  The 
factorization  problem  now  seems  to  be  that  of  factoring  the  Hermitieui 
matrices  (H*XH)  and  Z.  Equation  (27)  has  been  solved  for  certain 
special  cases. 


X.  h ^ (HF)jn(^n^v) 

i _ n _ 
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Case  1.  - In  the  first  case,  H is  a column  matrix  and  the  z*s 
are  independent.  This  is  the  case  solved  herein  (eq.  (12))  in  which 
only  the  corresponding  row  of  F (emd  C)  is  determined.  It  is  noted 
that,  in  general,  if  any  row  of  C vanishes,  then  the  corresponding 
row  of  F likewise  vanishes. 


Case  2.  - In  the  second  case,  H is  a diagonal  matrix  and  the  z’s 
are  independent . The  general  solution  for  this  case  is 


rv 


^r,l 

HrZy- 


(28) 


where 


Hj.r(zyWr) 

Ayj.  — — — 

H^Zy 

Frrl^  = Prf 

(ZyZy)  = 1^1^ 

The  functions  Hj.  and  Zy  are  minimum  phase  structures  and 

Ayj.  = Ayr,l  \r,2  where  Ayj.^3_  and  stable  structures. 

Case  5.  - Every  element  of  H and  is  stable  and  the  z’s 

are  independent  in  the  third  case.  The  solution  for  this  case  is 


rv 


-z 


(H  ^)rs'^sv,l 


where 


(wgZv) 


‘sv 


Try  is  defined  in  Case  2,  and  Agy  is  to  be  separated  into  the  two 
parts  Asv,l  + Asv,2  as  in  Case  2. 

Case  4.  - In  the  fourth  case,  = BDjMjjj  and  the  z's  are 
independent  where  HB  = DjDj  = 1 and  every  element  of  M and  M“^  is 
stable.  The  solution  for  this  case  is 


2789 


2789 
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Z(M~  )rs^sv,l 

Zv 


where 


BDs(wsZv) 


Here  is  defined  in  Case  2,  and  Ag^  is  to  be  separated  into  the 

two  parts  Ag^^^  + Ag^^2  Case  2. 
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Figure  1.  - General  linear  closed-loop  controlled  system;  one 
manipulated  variable.  Circle  indicates  operation  of  addition. 


I 1 


Figure  2.  - Representation  of  general  linear  controlled  system 
using  only  stable  boxes  E,  F,  and  G;  one  manipulated  vari- 


■ - E 

J 


able;  stable  if  G 
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Figure  3.  - Open-loop  system  equivalent  to  closed- 


loop  system  of  figure  2 vhen  G = 


■E 


Vj- 


Figure  4.  - Illustration  of  example  1,  speed  con- 
trol of  turtojet  engine. 
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Figure  5.  - Illustration  of  example  2,  spjeed  control  of  turbo- 
jet engine  vlth  temperature  variation  constrained. 


I 


Figure  6.  - Frequency  response  of  engine  speed  to  fuel  pressure  of  axial -flow 
turbojet  engine. 
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Figvire  7.  - Integral-square  errors  of  controlled  axial-flow  turbojet  engine  for 
range  of  loop  gains  and  Integral  time  constants. 


Figure  8.  - General  linear  closed -loop  controlled  system. 
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